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ABSTRACT 

^Steady  two-dimensional  flow  past  a  sail  is  considered.  The  sail  is 
assumed  to  be  supported  by  two  masts.  The  flow  and  the  shape  of  the  sail  are 
determined  as  functions  of  the  direction  o  of  the  flow  at  infinity  and  the 
Weber  number  The  full  nonlinear  problem  is  formulated  as  an  inteqro- 

differential  equation  for  the  shape  of  the  sail.  This  equation  is  discretized 
and  solved  numerically  by  Newton's  method.  Sail  profiles,  the  slack  in  the 
sail,  and  the  lift  coefficient  are  presented  for  various  values  of  o  and 
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SIGNIFICANCE  AND  EXPLANATION 


We  consider  the  flow  past  a  two-dimensional  sail  supported  by  two 
masts.  Previous  investiqators  (Thwaites1 ,  Nielsen2,  Tuck  and  Haaelgrove^, 
Vanden-Broeck  and  Keller^)  derived  approximate  solutions  for  this  problem  by 
linear izinq  the  boundary  condition  on  the  sail.  In  this  paper  we  provide  an 
exact  numerical  solution  of  the  fully  nonlinear  problem.  Ibis  solution  can  be 
used  to  estimate  the  range  of  validity  of  the  previous  approximate  solutions. 

Although  this  problem  is  very  specific,  the  technique  used  to  solve  it 
can  be  viewed  as  a  general  method  of  the  theory  of  aerodynamic  surfaces  formed 
by  flexible  membranes.  Examples  of  such  surfaces  are  parawings,  rotor  blades 
and  glider  winqs  utilizing  flexible  lifting  surfaces.  Other  applications 
involve  the  study  of  pneumatic  structures  such  as  those  sometimes  used  to 
cover  sport  arenas. 

Our  method  to  treat  the  basic  nonlinear  equation  is  new.  In  addition  it 
is  easier  to  use  than  the  approximate  methods  of  Thwaites1  and  Nielsen2,  which 
the  present  work  improves. 


The  responsibility  for  the  wordinq  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


NONLINEAR  TWO-DIMENSIONAL  SAIL  THEORY 


Jean-Marc  Vanden-Broeck 

1 «  Introduction 

We  consider  the  deformation  of  a  two-dimensional  sail  due  to  the 
steady  potential  flow  of  an  incompressible  inviscid  fluid  around  it  (see 
Fig.  1).  The  sail  is  supported  by  two  masts  and  is  characterized  by  its 
constant  tension  a  and  by  the  distance  c  between  the  two  masts.  The 
fluid  has  density  p  and  velocity  U  at  infinity.  As  we  shall  see,  the 
shape  of  the  sail  is  determined  hy  the  direction  ct  of  the  flow  at 
infinity  and  by  the  Weber  number 

X  •  2 pctl2/o  .  (1) 

1  2 

This  problem  was  first  considered  by  Thwaites  and  Nielsen  .  These 

authors  obtained  approximate  solutions  for  a  small  by  using  thin  aerofoil 

3 

theory.  Their  work  was  further  generalized  by  Tuck  and  Haselgrove  and 
Dugan*.  More  recently  VSnden-Broeck  and  Keller’’  obtained  an  asymptotic 
solution  for  X  small  for  arbitrary  values  of  a. 

In  the  present  paper  we  solve  the  fully  nonlinear  problem 
numerically.  In  Sec.  II  we  formulate  the  problem  as  an  integrodifferential 
equation  for  the  unknown  shape  of  the  sail.  In  Sec.  Ill  we  present  a 
method  to  solve  this  equation  numerically.  The  method  involves 
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discretization,  Which  converts  the  equation  into  a  set  of  nonlinear 
alqehraic  equations.  Then  it  employed  Newton's  method  to  solve  these 
equations.  The  results  obtained  are  discussed  in  Sec.  IV. 

Many  different  families  of  solutions  exist.  Tuck  and  Haselgrove 's6 
stability  analysis  shows  that  only  one  family  of  solutions  is  stable  for 
o  small.  Therefore  we  present  only  numerical  results  for  this  family.  It 
is  found  that  for  each  value  of  a  there  is  a  maximum  value  of  the  Weber 
number  X  above  which  this  family  fails  to  exist. 


II .  Formulation. 

Let  us  consider  the  steady  two-dimensional  flow  of  an  inviscid  in¬ 
compressible  fluid  past  an  inextensible  and  flexible  sail.  The  flow 
configuration,  the  sail  and  the  coordinates  are  sketched  in  Fig.  1. 

The  leading  edge  or  luff  of  the  sail  is  attached  to  a  mast  at  x  =  0, 
v  =  0  and  the  trailing  edge  or  leach  is  attached  to  another  mast  at 
x  =  c,  y  =  0.  At  infinity  we  require  that  the  flow  have  speed  U  and 
direction  a. 

We  introduce  dimensionless  variables  by  taking  c  as  the  unit  length 
and  U  as  the  unit  velocity.  Let  the  shape  of  the  sail  be  described  in 

dimensionless  variables  by  y  *  F(x),  0  <_  x  <_  1.  The  conditions  of 
attachment  to  the  masts  imply 

F (0)  =  F(l)  =  0  .  (2) 
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The  difference  between  the  pressures  on  the  two  sides  of  the  sail  is 


balanced  by  ie(x)a,  where  k  is  the  curvature  of  the  sail  and  a  is  the 
constant  tension  in  it.  By  using  Bernoulli's  equation  we  obtain  in  di¬ 
mensionless  variables 

<(x)  =  j  [q'(x)2  -  q+(x)2]  .  (3) 

Here  q+(x)  and  q  (x)  are  the  flow  velocities  above  and  below  the  sail, 
respectively.  The  p irameter  A  is  the  Weber  number  defined  by  (1). 

It  is  convenient  to  introduce  the  complex  perturbation  velocity 

a  -  ib  =  u  -  cos  a  -  i(v  -  sin  a)  .  (4) 

Here  u  and  v  denote  respectively  the  chordwise  and  normal  components 
of  the  velocity.  The  function  a  -  ib  is  an  analytic  function  of 
z  =  x  +  iy  which  vanishes  at  infinity.  Therefore  we  can  write,  by  using 
the  Plemelj  formula 

a+(x)  +  a~(x)  -  i[b+(x)  +  b  (x) ]  ^ 

)  -  a+(xx)]  -  i[b'(Xl)  -  b+(Xl)]}[l  ♦  iF ' (Xj ) ]  ^ 

(Xj  -  x)  +  i [F (Xj  -  F(x)]  dXl 

0 

The  integral  in  (5)  is  of  Cauchy  principal  value  form. 

The  velocity  must  be  tangent  to  the  sail  on  both  sides.  Thus  we 

have 

v4(x)  *  u4(x)F'(x)  ,  0  x  <  1  .  (6) 

Using  the  definition  (4),  we  can  rewrite  (6)  in  the  more  convenient  form 
b+(x)  -  b"(x)  =  [a+(x)  -  a"(x)]F'(x)  ,  (7) 


.  i  I  f[a'(xi 

TT  ! 
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b+(x)  +  b~(x)  *  -2  sina  ♦  [a+(x)  +  a‘(x)]F'(x)  ♦  2cosa  F'(x)  .  (8) 

By  substituting  (7)  into  (5)  and  taking  real  and  imaginary  parts  we 
obtain 


1  r  + 


b+(x)  ♦  b"(x)  =  ± 


[a  (Xj)  -  a  (x1)][l  +  F'CXj)  ] (Xj  -  x) 
(Xj  -  X)2  *  [F(Xj)  -  F(X)]2 


dxj  ,  (9) 


+  1  r  [a+(Xj)  -  a  (x  )][1  +  F'CXj)2] [FfXj)  -  F(x)] 

a  (x)  ♦  a  (x)  ■  -  —  I  - 5 - 2 - dxl 

71  J0  (Xj  -  x)*  *  [F(Xj)  -  F(x)] 


(10) 


We  now  rewrite  (3)  in  the  form 


F"(x) 


[1  +F.<x>V'2  4 


+  4  Ia+(x)  -  a”(x) J [ 1  +  F'(x)2][2  cob  a 


(11) 


+  a  (x)  +  a  (x )  1  ■  0  . 

Finally  we  impose  the  Kutta  condition  at  x  =  1  by  requiring 


a+(l)  =  a"  (1)  . 


(12) 


The  condition  (12)  ia  clearly  not  satisfactory  when 

F*  (x)  ~  ( 1  -  x)v  with  v  <  -  y  aa  x  ♦  1  . 

However  such  exceptional  cases  did  not  occur  in  the  numerical  solution.  He 
shall  consider  only  values  of  a  between  0  and  V 2  so  that  the  Kutta 
condition  is  imposed  at  the  trailing  edge. 

For  given  values  of  A.  and  a  relations  (8) -(12)  define  a  nonlinear 
system  of  integrodifferential  equations  for  the  unknowns  F(x),  a+(x), 
a~(x),  b+(x)  and  b~(x). 
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By  substituting  (7),  (9)  and  (10)  into  (8)  and  (11)  we  obtain  the 
following  reduced  system  of  equations  for  the  two  unknowns  a+(x)  -  a  (x) 
and  F(x): 

j  [1  +  F'(x1)2][a+(x1)  -  a”(x1)]{(x1  -  x)  ♦  F'(x)[F(x1)  -  F(x)]} 
j,  (Xj  -  X)2  +  tF(Xj)  -  F(X)]2 

=  -2sim  +  2cosa  F*(x)  ,  (13) 


- "”~Xy~7T2  ♦  j  ta+(x)  -  a"(x)][l  +  F'(x)2] 

[1  ♦  F’(x)Y/Z 


i 


1  r  + 


2cosa  -  — 

IT 


[a+(Xl)  -  a~  (xx)]  (1  ♦  F'(x1)2][F(x1)  -  F(x)] 
(Xj  -  x)2  +  [F (Xj)  -  F(x)] 2 


dxi| 


=  0  . 


(14) 


This  system  is  solved  numerically  in  the  next  section. 

From  the  solution  we  can  then  compute  the  slack  l  in  the  sail  by  the 
formula 


4  = 


(1  +  [F’(x)]2)1/2dx  -  1  . 


(15) 


Following  Thwaites1  we  define  the  parameter  B  by  the  relation 


1/2 

B  =  sin2a/24  ' 


(16) 


The  lift  coefficients  is  obtained  from  the  solution  by  substituting  (7) 
into  the  formula  (15)  given  by  Vanden-Broeck  and  Keller^.  Thus  we  obtain 


CL  = 


[a  (x)  -  a+(x)]{l  +  [F'(x)]2}dx 


(17) 


-5- 


III.  Numerical  analysis. 

In  order  to  solve  the  system  of  equations  defined  by  (2),  (12),  (13) 
and  (14)  we  introduce  the  mesh  points 


Xj  =  (I  -  1)/(N  -  1)  , 

I  =  1, . . . ,N  , 

(18) 

and  the  corresponding  unknowns 

FI  =  F(xi)  * 

I  =  1,...,N  , 

(19) 

Ai  =  a+(xP  '  a"(xj)  » 

I  =  2.....N  . 

(20) 

Relations  (2)  and  (12)  imply  =  0  and  *  0  so  that  there  are 

only  2N  -  4  unknowns  Fj  and  Aj. 

We  shall  satisfy  the  equations  (13)  and  (14)  at  the  N-3  intermediate 
mesh  points 

xI+1/2  =  (xI  +  XI+1)/2  *  I  -  2.....N-2  .  (21) 

First  we  compute  F'(Xj),  F'(xI+1/2),  F"(xI+1/2),  F(xI+1/2)  and  a+(xI+1/2) 
-  a  (Xj+1/2)  *n  terms  °F  Fj  and  Fj  by  four  point  formulas. 

Next  we  evaluate  the  integrals  in  (13)  and  (14)  by  using  the 
trapezoidal  rule  with  the  mesh  points  Xj.  The  mesh  points  Xj  are  locally 
symmetric  about  XI+1^2  and  the  quadrature  formula  is  also  symmetric. 
Therefore  the  contribution  from  the  neighborhood  of  the  singularity 


cancels  out,  permitting  us  to  evaluate  the  Cauchy  principal  value  as  if 
it  were  an  ordinary  integral. 


In  doing  so  it  is  important  to  notice  that  the  integrands  in  (13) 

-1/2 

and  (14)  behave  like  x  as  x  -*■  0.  We  handle  these  singularities  by 
subtracting  the  singular  part  before  integrating  numerically.  We  then 
integrate  the  singular  part  analytically. 

By  discretizing  (13)  and  (14)  we  obtain  2N  -  6  nonlinear  algebraic 
equations  for  the  2N  -  4  unknowns  Fj  and  A”.  An  extra  equation  is 
obtained  by  imposing  the  Kutta  condition  at  x  =  1  via  the  relation 

lim  (1  -  x)1/2[a+(x)  -  a"(x)]  =  0  .  (22) 

x+1 

A  three  point  Lagrange  extrapolation  formula  is  used  to  evaluate  the 
limit  in  terms  of  A^ .  The  last  equation  is  obtained  by  imposing  (2)  by 
means  of  a  three  point  Lagrange  formula. 

For  given  values  of  X  and  a  this  system  of  equations  is  solved  by 
Newton's  method.  The  values  of  g  and  can  then  be  computed  by  quadra¬ 
ture  from  (16)  and  (17). 

IV.  Discussion  of  the  results. 

We  have  used  the  method  described  in  Sec.  Ill  to  compute  solutions  for 
various  values  of  X  and  a.  Many  different  families  of  solutions  exist. 
They  have  been  studied  for  a  small  by  Thwaites1.  These  families  can  be 
computed  for  arbitrary  a  and  X  in  the  following  way.  For  a  given  value 
of  X  we  compute  the  solution  for  a  small  value  of  a  by  using  the 
uniform  stream  (i.e.  the  solution  corresponding  to  a  =  0)  as  the  initial 
guess.  Providing  X  is  not  close  to  any  of  the  critical  values  described 
by  ihwaites1,  the  numerical  scheme  converges  rapidly  to  Thwaites'  solution 
for  the  given  value  of  X.  Once  a  given  solution  has  been  obtained, 
however,  a  type  of  'boot-strap*  technigue  is  employed,  that  is  a  converged 
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solution  for  one  value  of  X  is  used  as  the  initial  quess  for  a  solution 
with  X  altered  by  a  few  per  cent.  Similarly,  for  fixed  X,  the  anqle 

a  could  be  varied  usinq  a  neiqhborinq  solution  as  a  startin'*  quess. 

Among  all  these  possible  solutions,  only  one  is  "reasonable"  in  the 

sense  that  it  has  a  simple  unimodal  shape  cambered  to  lee  of  the  wind . 

Tuck  and  Haselgrove* s^  stability  analysis  shows  that  this  family  of  solu¬ 
tions  is  the  only  stable  one  for  a  small.  Therefore  we  extend  only  this 
family  to  the  nonlinear  case. 

Typical  profiles  are  shown  in  Figs.  1  and  2.  The  values  of  3  and 
for  a  =  ir/4  are  shown  as  functions  of  X  in  Figs.  3  and  4.  The  broken 
lines  correspond  to  the  asymptotic  formulas  (14)  and  (20)  given  by 
Vanden-Broeck  and  Keller5.  The  asymptotic  formula  for  3  agrees  with  the 

numerical  results  within  fifteen  percent  for  X  <  0.3. 

In  Piq.  5  we  present  the  values  of  X  versus  3  for  a  -  */100  an<J 

a  *  w/10.  The  broken  line  corresponds  to  Thwaites*  solution.  It  shows 
that  his  solution  is  a  qood  approximation  for  a  small,  but  it  is  not 
uniform  as  B  tends  to  zero.  The  discrepancy  between  his  solution  and  the 
exact  numerical  solution  increases  rapidly  as  3  qets  close  to  zero.  In 
the  particular  case  3-1,  the  values  of  were  computed  for  various 

values  of  a.  The  numerical  values  of  CL  were  found  to  aqree  with 
Thwaites*  solution  within  5%  for  a  <  ir/18,  20%  for  °  <  n/6  and  30%  for 
a  <  */4. 

Pigs.  3  and  5  and  similar  results  obtained  for  different  values  of  a 
indicate  that  for  each  value  of  a  there  exists  a  maximum  value  X^(a)  of 
the  Weber  number  above  which  there  is  no  solution  of  the  kind  considered 
here.  For  a  qiven  value  of  a  there  are  two  solutions  for  each  value  of 
the  Weber  number  in  an  interval  below  the  maximum  value.  This  non- 
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uniqueness  is  likely  to  he  of  mathematical  interest  only.  The  reason  for 
this  is  that  the  branch  on  the  left  of  the  maximum  is  probably  unstable. 
After  all  it  corresponds  to  an  increase  of  the  tension  arisinq  from  an 
increase  in  the  slack.  This  suqqests  that  the  physically  realizable 
solutions  correspond  to  X  >  X^(ot). 

As  a  increases,  the  slope  of  the  profile  increases  at  the  leading 
edge.  This  was  found  to  limit  the  accuracy  of  the  numerical  scheme. 
Accurate  solutions  for  ot  >_  tt/3  and  X  ~  XQ(a)  could  not  be  computed  with 
N  35. 

Finally  it  is  worthwhile  mentioninq  that  the  present  model  is  not 

since  separation  is  then  likely  to  occur.  A  more 
larqe  was  proposed  by  Duqan^. 


physical  for  a  larqe 
realistic  model  for  a 
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Figure 


Captions  for  Figurea 


Figure  1 .  Computed  profile  of  a  sail  with  a  =  w/6,  X  a  0.82  and 
B  =  3.0.  The  unit  length  in  y  is  c. 

Figure  2.  Computed  profiles  of  a  sail  with  o  a  w/4  of 
X  *  0.5,  0.8  and  1.0.  The  corresponding  values  of  8  are 
respectively:  6,  2.8  and  1.6.  The  unit  length  in  which  both  x  and 
y  are  measured  is  c. 

Figure  3.  Values  of  the  Weber  number  X  as  a  function  of  3  for 

u  «  w/4.  The  dashed  curve  corresponds  to  Vanden-Broeck  and  Keller's 
asymptotic  formula. 

Figure  4.  Values  of  the  lift  coefficient  CL  as  a  function  of  8 
for  a  a  w/4.  The  dashed  curve  corresponds  to  Vanden-Broeck  and  Keller's 
asymptotic  formula. 

Figure  5.  Values  of  the  Weber  Number  X  as  a  function  of  3  for 

a  a  w/ioo  (curve  a)  and  a  *  w/io  (curve  b).  The  dashed  curve 
corresponds  to  Thwaites '  calculations. 
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Steady  two-dimensional  flow  past  a  sail  is  considered.  The  sail  is  assumed 
to  be  supported  by  two  masts.  The  flow  and  the  shape  of  the  sail  are  determined 
as  functions  of  the  direction  a  of  the  flow  at  infinity  and  the  Weber  number  X. 
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the  shape  of  the  sail.  This  equation  is  discretized  and  solved  numerically  by 
Newton's  method.  Sail  profiles,  the  slack  in  the  sail,  and  the  life  coefficient 
are  presented  for  various  values  of  a  and  X. 
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